Detection of CO and ethane in Comet 21P/Giacobini-Zinner: 
Evidence for variable chemistry in the outer solar nebula. 
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Abstract 

Ethane and carbon monoxide were detected in a short-period comet of probable Kuiper belt origin. 
Ethane is substantially less abundant compared with Hyakutake and Hale-Bopp, two comets from 
the giant - planets region of the solar nebula, suggesting a heliocentric gradient in ethane in pre- 
cometary ices. It is argued that processing by X-rays from the young sun may be responsible. 
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1. Introduction 

Knowledge of the composition and 
structure of cometary nuclei is central to 
understanding the evolution of material in our 
natal interstellar cloud, and in the pre- 
planetary nebula (cf. Mumma et al. 1993, 
Irvine et al. 1999). A key question is 
whether the composition of pre -cometary ices 
varied with distance from the young sun. 

Oort cloud comets likely formed in the 
giant - planets region of the nebula, where 
temperatures ranged from -200 K near 
Jupiter to -40 K near Neptune. Kuiper belt 
comets formed in-situ beyond Neptune where 
temperatures as low as -20 K may have 
prevailed. Chemical variations could result 
from a number of (poorly constrained) 
processes, such as partial vaporization of icy 
grain mantles during nebular infall, 
condensation fractionation of nebular gas, 
thermal processing of ices, hydrogenation of 
icy grain mantles, energetic processing by X- 
rays, and UV photochemistry. The efficacy 
of such processing probably varied greatly 
with heliocentric distance and/or with time, 
enforcing corresponding differences in the 
volatile fractions of comets. For example, 
comets formed in the Kuiper belt 
(21P/Giacobini-Zinner is probably one 
example) may have incorporated unusually 
large amounts of apolar volatiles such as 
nebular CO, neon, methane, and acetylene 
(and ethane, if it existed in the nebular gas), 
if temperatures there were as low as 20 K. 
Icy grain mantles may have been more 
efficiently hydrogenated in regions of high 
H-atom density (e.g. an X-ray Dissociation 
Region), converting condensed phase CO to 
H 2 CO and CH 3 OH, and C 2 H 2 to ethane. 
Thus, the native abundances of these volatile 
species comprise a key test of cometary 
origins, and of nebular processing. 

Cometary abundances do seem to reflect 
local conditions in the comet-forming 


regions. A dichotomy in carbon chemistry 
was found in a study of 85 comets (A'Heam 
et al. 1995), suggesting that origin-related 
chemical differences are general. About one- 
half of the Jupiter - family comets (mainly of 
Kuiper-belt origin) were depleted in C 2 and 
C 3 relative to CN, compared with “typical” 
comets (including those formed in the giant - 
planets region), and A’Heam et al. suggested 
this was evidence of different processing in 
the trans-Neptunian region. The volatile 
composition of the long-period comets 
Hyakutake and Hale-Bopp (probably formed 
in the giant - planets region) is similar to icy 
grain mantles in dense interstellar cloud cores 
(Mumma 1996, Mumma 1997, Crovisier 
1999), and neither is depleted in C 2 
(Schleicher et al. 1997, Hicks and Fink 
(1997). 

The meaning of the C 2 depletion 
observed for some Jupiter-family comets is 
uncertain because the origin of C 2 is poorly 
understood. Combi and Fink (1997) 
critically reviewed the possible production 
mechanisms and suggested CHON grains as 
the likely source, based on the strict R’ 2 
dependence of the scalelength for the C 2 
precursor. Possible gaseous precursors 
(acetylene and ethane) were first detected in 
comet Hyakutake (Brooke et al. 1996, 
Mumma et al. 1996), and both species were 
later detected in comet Hale-Bopp (Mumma et 
al. 1997, Weaver et al. 1999). Sorkhabi et 
al. (1997) suggested acetylene (C 2 H 2 ) as the 
source in Hyakutake based on photochemical 
considerations and laboratory measurements. 
Photochemical production of C 2 from 
ethane(C 2 H 6 ) is less well studied. Molecular 
structure considerations suggest that the 
quantum yield for C 2 should be small but this 
may be compensated by the higher 
production rate of ethane in these comets (~ 
0.6 %). In the absence of direct 

measurements of C 2 H 6 and C 2 H 2 , the 
depletion of C 2 in some Jupiter-family 
comets can not be uniquely interpreted. 
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Giacobini-Zinner (hereafter G-Z) is one 
of the comets “depleted” in C 2 . It is in a 
short period (6.61 yr.), prograde, low 

inclination (31.86°) orbit, and is thought to 
have originated in the trans-Neptunian 
(Kuiper belt) region. Its last favorable 
apparition was in 1985, when observations 
confirmed earlier findings that it is depleted in 
C 2 and C 3 by factors of about 5-10, even 
though CN/OH was “normal” (Cochran and 
Barker 1987, Schleicher et al. 1987, Konno 
and Wyckoff 1989, Beaver et al. 1990). The 
present apparition presented the first 
opportunity to directly test the abundance of 
potential C 2 precursors. 

We report herein the first detections of 
ethane and CO in a Jupiter-family comet. 
When compared with two comets from the 
giant - planets region (Hyakutake and Hale- 
Bopp), our production rates and abundance 
ratios for G-Z imply that ethane was depleted 
in pre-cometary ices from the trans- 
Neptunian region of the solar nebula. 

2. Observations and Results 

We obtained pre-perihelion observations 
of comet G-Z spanning UT 1998 October 2- 
10, using the Cryogenic Echelle Spectrometer 
(CSHELL; Greene et al. 1993) at the NASA 
Infrared Telescope Facility, Mauna Kea, HI. 
During this interval, the heliocentric distance 
decreased from R = 1.25 au to 1.19 au, the 
geocentric distance decreased from A = 1.08 
au to 1.03 au, and the geocentric velocity of 
the comet was roughly constant at 
approximately -9.5 km s' 1 . CSHELL 
features a 256 x 256 pixel InSb array, a 30 
arc -second long slit which we oriented east- 
west on the sky, and a single-pixel spectral 
dispersion of ~2.7 km s' 1 . The plate scale is 
0.2 arc-seconds per pixel. At each grating 
setting, cometary data were acquired using a 
sequence of four scans (ABBA), with an 
integration time of 5 minutes on source per 
sequence of ethane scans (two minutes for 


CO). Images of the comet at 2.2 [im were 
taken before and after each sequence of 
scans, to monitor and correct for cometary 
drift. For each grating setting, spectra of 
infrared standard stars (obtained through a 4 
arc-sec wide slit) were acquired for absolute 
flux calibration of the comet spectra. 

The data were processed using algorithms 
specifically tailored to our cometary 
observations (cf. Dello Russo et al. 1998). 
Initial data processing included flat-fielding, 
removal of high dark current pixels and 
cosmic ray hits, and registration of spectral 
frames such that the spectral dimension falls 
along a row and the spatial dimension is 
orthogonal to this. Atmospheric 

transmittance models were obtained using the 
Spectrum Synthesis Program (SSP) (Kunde 
and Maguire 1974) which accesses the 
HITRAN-1992 Molecular Data Base 
(Rothman et al. 1992). SSP models were 
used to assign wavelength scales to the 
spectra, and to establish absolute column 
burdens for relevant absorbing species in the 
terrestrial atmosphere. The transmittance 
model was binned to the instrumental 
sampling interval, convolved to the resolution 
of the comet spectrum and normalized to the 
cometary continuum. The extracted comet 
spectra reveal emission from two Q-branches 
of C 2 H 6 (Fig. 1A) and two emission lines 
from CO (Fig. IB), all at the correct Doppler- 
shifted frequencies. 

On October 9.2 and 10.3 we observed 
emission from the r Qo and PQi branches of 
the V7 band of ethane (Fig. 1A). The C 2 H 6 
observations employed a 2 arc-second slit 
width resulting in a spectral resolving power 
v / Av ~ 1.5 x 10 4 . The telescope was 
nodded 15 arc-seconds along slit (both A and 
B beams contained comet signal) to cancel 
background sky and telescope emission. On 
October 2.2, we detected emission from the 
R0 and R1 lines of the CO v=l-0 
fundamental band near 4.7 /im. We used the 
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1 arc-second wide slit, providing v / Av - 
2.1x1C)4. This higher resolving power was 
necessary for separating the cometary 
emissions from corresponding absorption by 
CO in the terrestrial atmosphere, given the 
relatively small geocentric Doppler shift (Fig. 
IB). For CO, the telescope was nodded 60 
arc-seconds between positions A (comet in 
slit) and B (blank sky), so as not to cancel 
possible extended CO emission. 

We calculated molecular production rates 
from the line emission contained within 
nucleus-centered spectral extracts (2x2 arc- 
second for C 2 H 6 , and 1x3 arc-second for 
CO), as for comet Hyakutake (cf. Mumma et 
al. 1996). However, our work with comet 
Hale-Bopp demonstrated that these nucleus- 
centered rates are smaller than the true 
production rates for native release by 
approximately a factor of two, largely owing 
slit losses deriving from the point-spread- 
function (cf. Dello Russo et al. 1998 for 
details). The production rates presented 
below were scaled upward by two to account 
for this effect. 

The extraction of total production rates 
requires knowledge of the rotational 
temperature when only a portion of the 
molecular band is measured. We detected too 
few lines in G-Z to permit determination of a 
self-consistent rotational temperature, so we 
adopt a value (35 K ± 10 K) that is typical of 
comets with similar production rates and 
observed under similar conditions. Bockelee- 
Morvan et al. (1994a) measured rotational 
temperatures for methanol in comets C/1989 
XI (Austin) (Q(H 2 0) ~ 4 x 10 28 s’ 1 , T rot = 

24 ± 8 K) and C/1990 K1 (Levy) (Q(H 2 0) ~ 

25 x 10 28 s' 1 , T rot = 23 ± 15 K, 41 ± 10 K, 
and 24 ±7 K). They recommend a 
temperature range 20 -40K for methanol in 
these two comets. Methanol rotational 
temperatures in 22P/Kopff (Q(H 2 0) ~ 4 x 
10 28 s' 1 , T rol = 27 ± 12 K) (Biver 1999), and 
1 09P/Swift-T uttle (Q(H 2 0) - 50 x 10 28 s' 1 , 


Trot = 45 ± 7 K) (Bockelee-Morvan et al. 
1994b) are consistent. The shape of the 3.52 
^m v 3 band of methanol in Swift-Tuttle was 
best reproduced by a rotational temperature of 
~ 50K (Davies et al. 1993). Rotational 
temperatures obtained for water in 1 P/Halley 
(Q(H 2 0) - 22 x 10 28 s' 1 , T rot -40 K and 
Q(H 2 0) - 140 x 10 28 s' 1 , T rot - 50 K) are 
also consistent (Mumma et al. 1986, Weaver 
et al. 1987, Bockelee-Morvan and Crovisier 
1987). 

Comet Hyakutake at closest approach 
(0.1 au) provided T rot -70 K for water, 
methanol, and formaldehyde, but Hyakutake 
is not a good analogue for G-Z since the 
measured temperatures refer to a coma region 
much closer to the nucleus owing to the small 
geocentric distance. Hale-Bopp, with its 
huge water production rate, is also not a good 
analogue for G-Z since thermalization of fast 
H-atoms is more efficient, leading to higher 
temperatures. Optical depth effects in the far 
infrared water cooling lines are also far 
stronger, again leading to higher temperatures 
(Bockelee-Morvan 1987, Bockelee-Morvan 
and Crovisier 1987). We assume that 
methanol, water, and ethane attain similar 
rotational temperatures, since electron 
collisions seem to control rotational excitation 
in the intermediate coma (Xie and Mumma 
1992, Bockelee-Morvan et al. 1994a, 
Bockelee-Morvan et al. 1994b). We 
therefore adopt the range 25 - 45 K as 
reasonable for both CO and ethane in G-Z, 
with a most likely value of about 35 K, and 
we evaluate production rates for 25, 35, 45, 
and 50 K (Table 1). 

Using a fluorescence efficiency (at 1 au 
and 35 K) of 10.72 xlO 5 photons s 1 for the 
combined r Qo and PQi branches and the 
measured combined line flux 1(1.0 ± 0.2) x 

— 1 8 9 L V / 

10 watts m'"], we find Qc 2 H6 = (3.5 ± 
0.8) x 10“ 5 molecules s 1 for the nucleus- 
centered extract. Including the correction 
factor of two, we obtain Qc 2 H6 = (7.0 ± 1.5) 
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x 10 25 molecules s' 1 . (Note: these two Q- 
branches were mislabeled in Mumma et al. 
1996, and inappropriate g-factors were used. 
However, two other errors were made that 
largely compensate this error. Revised 
production rates are in preparation, and the 
revised mixing ratio is given in the 
Discussion. The parameters used here are 
correct.) 

The line flux measured for the RO line of 
CO was (2.27 ± 0.56) x 10" 18 watts m’ 2 and 
the fluorescence efficiency (at 1 au and 35 K) 
is 1.53 x 10' 5 s’ 1 , leading to Q = (1.68 ± 
0.42) x 10 27 molecules s’ 1 after correcting 
for atmospheric transmission (0.70). The 
line flux measured for the R1 line was (1.97 
± 0.61) x 10' 18 watts m’ 2 and the 
fluorescence efficiency is 2.33 x 10’ 5 s’ 1 (at 1 
auand35 K), implying Q = (1.56 ± 0.48) x 
10 27 molecules s' 1 after transmittance 
correction (0.43). The weighted mean of the 
two lines provides Q C o = (3.28 ± 0.64) x 
10 27 molecules s' 1 , after including the 
scaling factor of two. The small aperture 
emphasizes the contribution from native CO, 
though a small contribution from an extended 
source is possible (cf. DiSanti et al. 1999a). 
For the rotational temperature of 35 K, the 
mixing ratio of native ethane to native CO is 
then (2.1 ± 0.6) x 10' 2 . The mixing ratio is 
not very sensitive to the assumed rotational 
temperature (Table 1). 

Cometary water is the paradigm with 
which all other native (parent) volatiles are 
compared. Its abundance can be estimated by 
indirect means (e.g. from radio and 
ultraviolet OH measurements, and from 
visible O^D) emission). Water can also be 
detected directly from the ground by infrared 
non-resonance fluorescence (Mumma et al. 
1995), but the need for frequent imaging 
limited our time on source and constrained 
application of this approach to G-Z. A line of 
H 2 0 (1)| — 1 io in the v 3 — v 2 hot-band, Dello 
Russo et al. 2000) should appear near the R1 


line of CO, at a Doppler-shifted frequency of 
2151.26 cm 1 . Although no strong line is 
seen, there is an excess flux at the proper 
frequency and an upper limit to the water 
production rate can be obtained. The 
measured line flux is (8.9 ± 4.8) x 10" 19 
watts m’ 2 (giine = 5.1 x 10' 7 s’ 1 at 1 au and 35 
K), corresponding to Qh 2 o = (3-2 ± 1.7) x 
10 28 molecules s' 1 (a transmittance correction 
of 0.87 and scaling by two are included). 
The water production rate on Oct. 2.2 was 
less than 6.6 x 10 28 molecules s' 1 , at the 
95% confidence limit. 

Water production rates were obtained by 
several indirect methods during the 1985 
apparition (McFadden et al 1987, Schleicher 
et al. 1987, Gerard et al. 1988, Konno and 
Wyckoff 1989, Tacconi-Garman et al. 1990, 
Combi and Feldman 1992). Measurements 
were taken at larger and smaller heliocentric 
distance, but not at the exact distance of our 
observations. Based on an interpolation of 
those measurements, the water production 
rate expected for G-Z during early October 
1998 is ~ 5 x 10 28 s' 1 . Our measured 
production rate is consistent with this 
estimate, at the one-sigma confidence level. 
However, little can be inferred from this 
agreement, considering the low S/N ratio of 
our water detection. Also, the behavior of 
the comet may have differed in 1985 and it 
was not sampled at our exact heliocentric 
distance. 

Comet Giacobini-Zinner was also 
observed by Weaver et al. (2000) using 
CSHELL at IRTF (October 25 - 29, 1998), 
several weeks after the observations reported 
here. They did not detect either ethane or 
CO, but they did detect water vapor. Their 
three-sigma upper limit for ethane represents 
a depletion by ~ 10-fold (relative to water), 
compared with comets Hyakutake and Hale- 
Bopp. Their three-sigma upper limit for CO 
(~ 2-3 % relative to water) is considerably 
smaller than the abundance retrieved by us 
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(~10 % relative to water). Their non- 
detections of CO and ethane are surprising in 
view of the clear detections shown in Fig. 1 . 

3. Discussion 

It is interesting to compare production 
rates for hypervolatiles in the Jupiter-family 
comet (G-Z) with those for comets from the 
giant - planets region. The preliminary 
production rates reported for Hale-Bopp 
(Mummaet al. 1997) were revised following 
more mature data reductions that account for 
slit losses and also separate the native and 
distributed sources for CO (ethane, Dello 
Russo et al. 1999a; CO, DiSanti et al. 
1999a). The relative production rates for 
native ethane and native CO are C 2 H 6 : CO = 
0.2 : 10 for G-Z, compared with 0.6 : 12 for 
Hale-Bopp. Their mixing ratio is (2.1 ±0.6) 
x 10' 2 in G-Z (at 35K), compared with (4.7 
± 0.4) x 10" 2 in Hale-Bopp. If the mixing 
ratios were constant in G-Z during the period 
Oct. 2-10, the abundance ratio of apolar 
volatiles (C 2 H 6 : CO) in these comets differs 
significantly — in G-Z being less than one - 
half that in Hale-Bopp. The production rate 
of ethane relative to water in this Jupiter- 
family comet is also low compared with both 
long-period comets (C 2 H 6 : H 2 O: = (2.2 ± 
1.3) x IO - 3 in G-Z, (5.6 ± 0.4) x 10* 3 in 
Hale-Bopp, and (6.4 ± 1.5) x 10' 3 in 
Hyakutake (Dello Russo et al. 1999b)). A 
comparison of native CO production among 
these three comets will be possible after our 
ongoing analysis of extended and distributed 
sources in Hyakutake is complete. 

What could cause this depletion of ethane 
relative to CO in G-Z? It seems unlikely that 
selective loss of ethane during the comet’s 
active phase could change the C 2 H 6 /CO ratio, 
since both species are apolar hypervolatiles 
and should be affected in similar fashion. G- 
Z was observed at 1.2 au heliocentric, well 
within the distance at which water 
sublimation controls cometary activity, so it 


seems unlikely that heliocentric fractionation 
effects could affect the relative production 
rates. A simple view is that the ices in G-Z 
are depleted in C 2 H 6 , and that this difference 
is cosmogonic. It follows that some process 
affected C 2 H 6 differently in the trans- 
Neptunian region compared with the giant - 
planets region of the nebula. 

Mumma et al. (1996) argued that the 
abundances in Hyakutake were consistent 
with conversion of condensed phase C 2 H 2 to 
C 2 H 6 and of CO to H 2 CO and CH 3 OH, 
through H-atom addition reactions on grain 
surfaces (cf. Tielens and Allamandola 1987, 
Chamley et al. 1995). The 5-10 fold 
depletion of C 2 in G-Z suggests a large 
current deficit in C 2 H 2 , if acetylene is its 
precursor. If C 2 H 6 in G-Z was produced by 
hydrogenation of C 2 H 2 , then its lower 
C 2 H 6 /CO ratio and low current C 2 H 2 
abundance could result from a lower initial 
abundance of acetylene in pre-cometary ice 
coupled with a high conversion efficiency of 
acetylene to ethane in the trans-Neptunian 
region. (The possible generation of C 2 from 
CHON grains in G-Z is probably negligible 
since the small (sub-^tm sized) grains 
responsible for the continuum at optical 
wavelengths are underrepresented. 
However, a firm conclusion regarding C 2 H 2 
must await its direct measurement on a future 
apparition.) 

The production of H-atoms in nebular gas 
and solids undoubtedly varied with distance 
from the young sun, and with time. It is 
known that Class I and II T-Tauri stars are 
very intense X-ray emitters (cf. Neuhauser 
1997), and it is likely that the young sun 
passed through a similar stage before entering 
the main sequence (Shu et al. 1997). The 
principal loss mechanism for these X-rays is 
photo-ionization of nebular H 2 , a process 
which ultimately produces H-atoms in great 
numbers through direct and dissociative 
recombination reactions (producing several 
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hundred H-atoms per KeV) (Glassgold et al. 
1997). Before nebular clearing in the giant - 
planets region, the increasing attenuation of 
this exciting X-ray flux with greater 
heliocentric distance may have caused a 
significant reduction in hydrogenation 
efficiency for ices in the trans-Neptunian 
nebular region. As nebular clearing 
proceeds. X-rays will penetrate more deeply 
and ionization fractions and H-atom densities 
will increase. Thus, the ratio of C 2 H 2 to 
C 2 H 6 in Kuiper belt comets could be larger 
or smaller than that in Oort cloud comets, 
depending on the time and place of their 
formation. Comet formation is generally 
thought to proceed more slowly at large 
heliocentric distance. 

Weaver et al. (2000) report a non- 
detection of both CO and ethane in comet G- 
Z, and their derived upper limits are 
significantly smaller than the production rates 
inferred from our detections. This 
discrepancy is most easily explained if G-Z is 
internally heterogeneous in its chemical 
composition. During its 1985 apparition, G- 
Z showed marked asymmetry in water 
production about perihelion, with a rapid 
reduction in Q(OH) occurring just prior to the 
heliocentric distance at which Weaver et al. 
observed in 1998 (R h = 1.11 - 1.08 au). It is 
possible that the vent responsible for gas 


release during our observations became 
increasingly inactive after R h ~ 1.2 au owing 
to seasonal effects, and that thereafter gas 
was released from a different vent. The 
second vent could be depleted in the 
hypervolatiles CO and ethane. If so, the two 
data sets could provide evidence supporting 
the long-sought chemical heterogeneity 
within a comet nucleus. Such heterogeneity 
could result if cometesimals were scattered 
radially in the solar nebula before being 
embodied in the final nucleus or were formed 
at times when physical conditions (e.g. local 
H-atom densities) differed greatly in the 
nebula. The depletion of hypervolatiles could 
also have a physical origin if the two vents 
experienced different thermal processing. 
These ideas can be more rigorously tested by 
measuring the compositions of other comets, 
and of Giacobini-Zinner upon its return to the 
inner solar system in 2005. 
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Table 1. Production rates for Parent Volatiles in 21P/Giacobini-Zinner 


Trot 

CO 

C2H6 

H 2 O 

K 

, n 27 -1 
10 s 

in 27 -1 
10 s 

, n 27 -1 
10 s 

25 

2.7 ± 0.5 

0.062 ± 0.013 

43 ±23 

35 

3.3 ± 0.6 

0.070 ± 0.015 

32 ± 17 

45 

3.8 ± 0.8 

0.077 ± 0.015 

29 ± 16 

50 

4.2 ± 0.8 

0.080 ± 0.017 

28 ± 15 


C 2 H 6 : CO 

Ratio 

0.023 ± 0.007 
0.021 ± 0.006 
0.020 ± 0.006 
0.019 ± 0.006 


li Spectral extracts of Comet Giacobini-Zinner, smoothed by a five-pixel wide rectangular 
filter. A. Two Q-branches from the V7 band of ethane on UT 1998 October 9.2 and 
10.3. B. Line emission in the CO v=l-0 band on UT 1998 October 2.2. Vertical bars 
below the line assignments represent the expected Doppler-shifted line positions. The 
dashed line is the synthesized atmospheric transmittance, convolved to the instrumental 
resolution (see text). 
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